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Abstract 



We study the effect of non- vanishing masses and mixings among neutrino 
flavours on the detection of neutrinos from stellar collapse by a water Cerenkov 
detector. We consider a realistic frame- work in which there are three neutrino 
flavours whose mass squared differences and mixings are constrained by the 
present understanding of solar and atmospheric neutrinos. We also include 
the effects of high dense matter within the supernova core. We find that the 
number of events due to the dominant process involving electron-antineutrinos 
may change dramatically for some allowed mixing parameters. Furthermore, 
contributions from charged-current scattering off oxygen atoms in the detec- 
tor can be considerably enhanced due to flavour mixing; such events have a 
distinct experimental signature since they are backward-peaked. 

PACS numbers: 14.60.Pq, 13.15+g, 97.60.Bw 
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I. INTRODUCTION 



The sighting of a supernova, SN1987a in the large Magellanic cloud (LMC) ffl, led to 
great excitement, since, for the first time, the neutrinos from stellar collapse were detected 
by earth-based detectors @,[| . Unlike electromagnetic radiation, which takes a long time to 
emerge from the collapsing core, neutrinos provide direct information about core collapse. 
The direct observation of neutrinos from SN1987a by the Kamiokande (KII) [Q and the 
1MB || detectors forms the beginning of a new phase in neutrino astrophysics with far- 
reaching implications for particle physics. Since then, the KII detector has been upgraded 
with tremendous improvement both in size and resolution, and many new detectors like 
SNO and Borexino will begin taking data soon. 

Immediately after SN1987a, several authors analysed PHTU[ the neutrino events recorded 
by the KII and the 1MB detectors. While the number of events were not statistically sig- 
nificant enough to obtain quantitative information on the neutrino spectrum, there was a 
qualitative agreement between the predictions from the core collapse mechanism and the 
observations. The present situation with improved neutrino detectors affords quantitative 
analysis of the neutrino events if a supernova collapse were to take place in the near future. 
While the observational scenario is positive, there has also been much progress in under- 
standing the properties of neutrinos, namely their masses and mixings, through the analyses 



of solar and atmospheric neutrino puzzles [11]. Both the solar and atmospheric neutrino 



deficits that are observed conclusively point to the requirement of (at least) three neutrino 
generations and mixing among them. The masses and mixings are then constrained by the 
observed deficit in the neutrino fluxes |12|,13|. 



In this paper, we analyse in detail the signatures of neutrinos from stellar collapse. 
The analysis is confined to Type II supernovae (which occur when the initial mass of the 
star is between 8-20 solar masses) since the neutrino emission from these are significant 
enough to make reasonable predictions. We consider the emission of all three types of 
neutrino (and antineutrino) flavours. There exist many models of stellar collapse. The 
present understanding of neutrino emission involves dividing the neutrino emission into two 
distinct phases — the neutronisation burst and thermal neutrino emission. The number of 
neutrinos emitted during the burst phase is only a few percent of the neutrinos emitted 
during the thermal or cooling phase of the proto-neutron star. While only v e is emitted 
during the earlier phase, neutrinos and antineutrinos of all types are emitted during the 
final phase. Much of the binding energy of the neutron star is radiated away as neutrinos 
while a small fraction (few percents or less) is deposited in the shock wave that blows away 



the mantle flbj]. Detailed predictions for the luminosity and average energy as a function of 
time are available |T5| for neutrinos emitted during the burst and the cooling phase. We use 
these predictions as an input in our analysis. 

An important fact to note is that the neutrinos, which are produced in the high dense 
region of the core, interact with matter before emerging from the supernova. The presence 
of non-zero masses and mixing in vacuum among various neutrino flavours results in strong 
matter dependent effects, including conversion from one flavour to another. Hence, the 
observed neutrino flux in the detectors may be dramatically different for certain neutrino 
flavours, for certain values of mixing parameters, due to neutrino oscillations. The effect of 
mixing on the neutrino signal from supernovae has been analysed in detail before by several 
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authors [|T6HT8|] . The effect of masses and mixing on time-of-ffight information has been 



discussed in [ 19jp0 



Our analysis is based on the pioneering work of Kuo and Pantaleone [IB| where they 
include mixing among all three neutrino flavours. However, unlike all the previous analyses, 
we take into account the constraints on the neutrino mixing and masses imposed [p~2|JT3[1 
by solutions consistent with the solar and atmospheric neutrino puzzles. There are several 
possible solutions here, including a purely vacuum solution for the solar neutrino problem. 
We choose some typical (allowed) values for the mixing angles, to illustrate the possible 
effects of mixing. This is important from the point of view of integrating known constraints 
on the neutrino masses and mixings in order to obtain a realistic picture of neutrino emissions 
from supernovae. We analyse in detail the dependence of the recoil energy spectra on the 
mixing parameters, at water Cerenkov detectors. In Sec. II we review the theoretical 
framework of our calculation, including a reanalysis of the matter effects on the neutrino 
spectrum in the hot dense core. Sec. Ill highlights the various inputs — neutrino fluxes 
and cross-sections in the detector — that we have used in order to compute the event rate 
expected at water Cerenkov neutrino detectors. The numerical computation of the total 
number of events for the time integrated neutrino spectrum, and the effects of neutrino 
oscillation, are discussed in Sec. IV, where the results are discussed and summarised. Some 
well-known results, adapted to the present situation, are discussed in Appendices A and B 
for completeness. 



II. THREE FLAVOUR OSCILLATIONS IN HIGH DENSE MATTER 

In this section, we discuss the mixing among three flavours of neutrinos (or antineutrinos) 
and compute the electron (or antielectron) neutrino survival probability, P ee (or P ee ). We 
will explicitly show that these are the only relevant probabilities. While the theoretical 
details of this mixing are well-known, the effects of super-dense matter, such as is found in 
the stellar cores, are non-trivial. We shall also use this analysis to set our notation. 

The three flavour eigenstates are related to the three mass eigenstates in vacuum through 
a unitary transformation, 











= u v 








. U 3 . 



(1) 



where the superscript v on the r.h.s. stands for vacuum. The 3x3 unitary matrix, U v , 
can be parametrised by three Euler angles (u, <f>, if>) and a phase. The form of the unitary 
matrix can therefore be written, in general, as 

U v = U 23 (if;) x U phase x C/ 13 (0) x U 12 (u) , 

where is the mixing matrix between the i th and j th mass eigenstates with the mixing 

angle It has been shown that the expression for electron neutrino survival probability, 
integrated over the time of emission and of absorption, is independent of the phase and the 

They can be set to zero without loss of generality and we have 



third Euler angle if) pl|j22|1 
the following form for U v : 



3 



c<j> 



(2) 



where s$ = sin0 and = cos<ft, etc. The angles u and <fi can take values between and 
7r/2. Recently, the chooz collaboration set a laboratory limit on V e oscillations |27| that 
resulted in a strong limit, < 12°, on the (13) mixing angle ||28|| . However, a combination of 
solar and atmospheric neutrino data allows for both large and small angle solutions for the 
(12) mixing angle, u. The angle ip is large, typically of the order of ip ~ 45° (although this 
is not relevant here). These constraints will be imposed later on in our numerical analysis. 

The masses of the eigenstates in vacuum are taken to be /ii, /12 and /X3. In the mass 
eigenbasis, the (mass) 2 matrix is diagonal: 



M 2 




tiii[ + 



H\I + AM 2 



(3) 



where the mass squared differences are given by <5 2 i = ^2 — An an d £31 = ^3 — Ah- Without 
loss of generality, we can take <5 2 i and S31 to be greater than zero; this defines the standard 
hierarchy of masses. Neutrino oscillation amplitudes are independent of the first term so 
we drop it from further calculation. In the flavour basis, therefore, the relevant part of the 
mass squared matrix has the form, 



AM 2 = XT AM 2 U v] 

= 8 31 M 3l + 8 21 M 21 , 



(4) 



where 



M, 



31 



M, 



21 
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\ ~ c <j) s 4> s uj ~ s 4> S u c u) s <f> s uj 



The relevant matter effects may be included by a modified mass squared matrix, 



where 



AM^ = 5 3 iM 31 + S 21 M 21 + A(r)M A , 



M A 




(5) 



(6) 



(7) 
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and A{r) is given by 

A(r) = V2G F N e (r) x 2E , (8) 

which is proportional to the electron number density, N e (r), in the supernova core. Here r 
is the radial distance from the centre of the star. The detailed modifications due to matter 
effects are discussed in Appendix A. 

The maximum value of A occurs at the core and is approximately 2 x lO 7 ^ eV 2 , where 
E is the neutrino energy in MeV. The modification due to the matter dependence is similar 
to the case of solar neutrinos, although, unlike in the case of solar neutrinos, all flavours are 
produced in the supernova core. 

It is clear that the mass squared matrix is no longer diagonal in the presence of matter; we 
therefore diagonalise AM^ in order to determine the matter corrected eigenstates. This is a 
difficult problem in general for arbitrary values of £31 and £21 ■ These are however constrained 
by the limits on them given by the simultaneous analysis of solar and atmospheric neutrino 
problems, namely 1CT 3 < < 1CT 2 eV 2 and 821 < 1CT 4 eV 2 , so that 831 ~ <5 32 ; the value 
of A for energetic neutrinos (of a few MeV to tens of MeV) in the core is therefore several 
orders of magnitude greater than these mass-squared differences. The eigenvalue problem 
may thus be solved perturbatively, with the following hierarchy: v4(core) ^> £31 ^> 821- As 
a result, the electron neutrino undergoes two well-separated resonances when the value of 
A(r) approaches the two mass squared differences. Following Kuo and Pantaleone |16[], and 



using the above mass hierarchy, the matter mixing angle m is given by (see Appendix A 
for more details), 

+ o- £31 sin 20 

tan 2<f) m = - — . 9 

d 3 i cos 20 — A 

At the point of production inside the core, A (core) ^> 831, thus, m — > tt/2. This makes 
further calculations extremely simple, since the electron neutrino is produced as a pure 
^3) mass eigenstate in the core of the supernova. The survival probability of the electron 
neutrino is simply given by the projection of the ^3) mass eigenstate on to the |z/ e ) flavour 
state in the detector, after correcting for the Landau-Zener jumps which may occur in the 
stellar matter during propagation. The average survival probability of the electron neutrino 
is therefore given by 



3 

y \u v -\ 2 



jjm 

ej 



2 



sin 2 P 3 + cos 2 sin 2 u P 2 + cos 2 cos 2 u> Pi . (10) 



Here and u are the vacuum mixing angles defined earlier and Pi denote the Landau-Zener 
jump probabilities among the mass eigenstates, 

Pi = PhPi , (11) 
P 2 = P h (l - P t ) , (12) 

P3 = (1 - P h ) , (13) 

where Ph and Pi denote the jump probabilities at the higher and at the lower resonances. 
In Appendix B we show that for the parameters values relevant in the case of neutrinos 
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produced in the supernova core, Ph is actually very close to zero. Therefore to a good 
approximation, we may write, 



Pee = sin 2 0. (14) 

This result implies that the propagation of neutrinos is adiabatic in the high dense core. 
We will support this conclusion in Appendix B. The information given above is not enough 
to obtain completely the survival and oscillation probabilities of the individual flavours. 
However, since the detectors we are interested in do not separately detect and u T , this is 
sufficient for our analysis. We shall therefore use this form for the survival probability for 
the numerical results calculated in the next section. 

We now consider the case of V e propagation in high dense matter. The only change in 
this case is that the matter dependent term in the relevant part of the mass squared matrix 
has the opposite sign (to that in Eq. (0)), that is, 



A(r) = -V2 G F N e (r) x 2E . (15) 

The analysis goes through as in the case of v e propagation through matter and the mixing 
angle, 0, for antineutrinos in matter is given by, 

, , 5 31 sin20 

tan20 m = ttt—J - ( 16 ) 

d 3 i cos 2(j) + A 

On using the fact that A(core) ^> 631, we obtain (fi m — > in contrast to the solution 
4>m — > tt/2 for electron neutrinos. Thus V e is produced in the mass eigenstate, |Fi), in the 
core of the supernova. There are no Landau-Zener jumps to consider in this case since the 
resonance conditions are never satisfied unless the mass hierarchy is altered. The propagation 
is therefore adiabatic and the survival probability is obtained by simply projecting the \Vi) 
eigenstate on to the flavour eigenstate in vacuum (at the detector). The antineutrino survival 
probability is therefore given by, 

P ee = cos 2 0cos 2 uj , (17) 
where <fi and uj are as usual the vacuum mixing angles defined earlier. 



III. NEUTRINO FLUXES AND CROSS-SECTIONS 

We need the inputs of neutrino flux emission at the supernova, and neutrino cross-section 
at the detector, in order to obtain the event rates. We begin with a discussion of the neutrino 
fluxes. 



A. Neutrino fluxes 

Following Kuo and Pantaleone |H| , we denote the flux of various flavours of neutrinos 
and antineutrinos produced in the core of the supernova by Ff , where i denotes all the 
flavours. In particular we use the generic label for flavours other than v e and V e since 
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F o = F o =F o =F o =F o / lg N 

All these flavours are produced via the neutral-current (NC) pair production processes and 
therefore have the same flux for all practical purposes. However, the v e and V e fluxes are 
different from each other and the rest since they are produced not only by pair production 
but also derive contribution from charged- current (CC) processes. 

In the presence of matter, the flux emerging from the core undergoes changes due to 
oscillations as was discussed in the previous section. The flux reaching the detector from a 
supernova at a distance d from earth is reduced by an overall geometric factor of l/(47r<i 2 ). 
Apart from this, there is a further modification of the observed flux due to oscillations. The 
flux on earth, in the various flavours, is given in terms of the flux of neutrinos produced in 
the core of the supernova by, 

F = P F° + P F° + P F° 

= i£-(l-P«)«-i2) , (19) 

where we have made use of the constraint J2j Pij — 1 an d P e x denotes the probability of a 
flavour /i or r neutrino emerging as an electron neutrino. Since v^- and z/ r -induced events 
cannot be separated in water Cerenkov detectors, their combined flux on earth may be 
written as 



IF = F 4- F 

= 2F° + (l-P ee )(F° -F°) . (20) 

Note that flavour mixing does not affect the total flux. 

Similar expressions hold for antineutrino flavours with appropriate changes, that is, 

Fv e = Fl-(l-P ee )(Fl-F°); (21) 

and 

2F„ = 2F° + (1 - P ee ) {Fl - F°) . (22) 

Since P ee ^ P ee , in general, the mixing breaks the equality of the u x and T7 X fluxes. 

We use the luminosity and average energy distributions (as functions of time) as given 
in Totani et al. |l5j , based on the numerical modelling of Mayle, Wilson and Schramm [14 



The neutrino number flux is described, in a given time interval, At, as a thermal Fermi 
Dirac distribution, 

dE iV V/ (exp^/T,-) + 1) ' 1 } 

for neutrinos of flavour j and energy E at a time t after the core bounce. Here % refers 
to the time-bin, t — to + iAt. Hereafter, we set the time of bounce, to = 0. The overall 
normalisation, iVo, is fixed by requiring that the total energy emitted per unit time equals 
the luminosity, £j, in that time interval. 

The thermal distribution that we have used (where the chemical potential has been set 
to zero) results in a flux that shows a slower fall with energy, E, than the results of the 
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corresponding numerical model used in Ref . |[15j ; however, the effect of this on the event rates 
is small, of the order of a few percent. Typically, this distribution corresponds to an average 
neutrino energy or temperature ((Ej) ~ 3.151}) of (E(v e )) ~ 12 MeV, {E(y e )) ~ 16 MeV, 
and (E(u x )) ~ 24 MeV. Beyond about 1 second after the bounce, the average energies remain 
constant over the emission times of the supernova. However, the luminosities decrease, with 
very little emission beyond 10 seconds. Hence, in order to compute the event rates, we 
consider neutrino emission up to 10 seconds after bounce. The total emitted energy in all 
flavours of neutrinos up to this time is about 2.7 x 10 53 ergs, which is more or less equally 
distributed in all flavours. The number of neutrinos emitted in each flavour, however, is not 
the same since their average energies are different. 



B. Interaction at the detector 

The basic quantity we are interested in is the distribution of events in the detector as a 
function of the energy of the detected particle. In the case of a water Cerenkov detector, this 
corresponds to the detection of a charged lepton in the final state. Here we are concerned 
with detection of electrons (or positrons) with energy, E e . The various processes of interest 
therefore are the interactions of the neutrinos 

1. with electrons in water as targets: 

u e (v e ) + e~ -> ut(v e ) + e~ , £ = e, fx, r ; (24) 

2. with free protons in water as targets: 

I7 e + p -> e + + n ; (25) 

3. with oxygen nuclei in water as targets: 

v e + 16 -> e~ + 16 F , 

V e + 16 -> e + + 16 iV . (26) 



The cross-sections, da/dE e , for all these processes, except the ones on oxygen, are well 
known [ 23], 24"1 . The oxygen cross- sect ions have been taken from Fig. 1 of Haxton [25]. As 



the interactions on protons and oxygen nuclei are purely CC interactions, they involve only 
u e and V e . Reaction (1) involves both CC and NC interactions for v e and V e and only NC 
interactions for all other flavours. 

The 77 e p cross-section is the largest, being proportional to the square of the antineutrino 
energy. In terms of total number of events, therefore, water Cerenkov detectors are mostly 
dominated by V e events. However, the different interactions in the detector have distinct 
angular signatures; this may be used to distinguish them. The elastic electron cross-sections 
are forward peaked, especially for neutrinos with energies ^ 10 MeV , while the proton 

cross-section is isotropic in the lab frame. Finally, the CC v e (y e ) cross-section on oxygen, 
although having a rather large threshold of 15.4 MeV (11.4 MeV) [J25|, increases rapidly with 
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incoming neutrino energy and is somewhat backward peaked. The higher the temperature 
at which the neutrino is emitted, the larger is this backward peak; hence it may be possible 
to distinguish this contribution from the rest by angular resolution as well, especially if there 
is substantial mixing between u e (or V e ) and v x since the latter have a considerably hotter 
spectrum. 



The time integrated event rate, from neutrinos of flavour j and energy E, as a function 
of the recoil electron (or positron) energy is given by, 



where the flux distribution, dF(j)/dE includes the effects of mixing in the hot dense core 
and the index t refers to any of the various processes through which the neutrino j can 
interact with the detector. Here N t refers to the number of scattering targets (of either e, p 
or 16 0) that are available in the detector. The total number of events from a given flavour 
of neutrino in a given bin, k, of electron energy (which we choose to be of width 1 MeV) 
then is the sum over all possible processes integrated over the bin width of the event rate: 



In the next section, we shall use this formula to compute the time integrated event rates 
for neutrino scattering with and without mixing, in water Cerenkov detectors, as a function 
of the detected electron (or positron) energy, in order to examine the effects of neutrino 
oscillations on supernova neutrino fluxes. 



We compute the time integrated event rate at a prototype 1 Kton water Cerenkov de- 
tector from neutrinos emitted by a supernova exploding 10 KPc away. Results for any other 
supernova explosion may be obtained by scaling the event rate by the appropriate distance 
to the supernova and the size of the detector, as shown in Eq. (p7|) . We assume the efficiency 
and resolution of such a detector to be perfect. Including these effects does not change the 
results by more than a few percents, as we will see. In fact, the maximum variation is at low 
energies, close to the threshold, where the low detector efficiency leads to lower detection 
rates. 

The following constraints derived from solar and atmospheric neutrino observations are 
imposed. We begin with the constraints in the neutrino sector. Here, the angle u does not 
play a role. As stated earlier, the (13) mixing angle is severely restricted: <ft < 12°. The 
solar and atmospheric neutrino problems allow for a wider choice in 0. This restriction on 
comes mainly from the chooz experiment ||27|| . Since the survival probability, P ee depends 
only on this angle (P ee = sin 2 0); this implies that P ee < 0.05 and is thus very small. The 



C. Event rates 




(27) 




(28) 



IV. RESULTS AND DISCUSSION 
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observed dynamics of electron-type neutrinos, therefore, is completely driven by mu- and 
tau-type neutrinos produced in the supernova: 



o 



2i^(F° + F°). (29) 

Hence, due to mixing, the original electron neutrino flux is virtually replaced by the \i or 
r neutrino flux. The cross-sections at the detector increase with energy. Since the average 
energy of v x is of the order of ~ 24 MeV while that of v e is ~ 11 MeV, the effect of mixing 
and matter in the dense core is to dramatically increase the number of events due to e-type 
neutrinos while reducing the corresponding v x contribution. 

We now discuss the antineutrino sector. While the same limits apply on <p, we now have 
to consider the limits on the (12) mixing angle to as well. The constraints on uj mainly emerge 



from the solar neutrino problem. (For a recent review, see Ref. [30]). The best global msw 
fit gives 5 \2 ~ 10 -6 eV 2 and sin 2 2uj = 5.5 x 10~ 3 . There is also a large angle solution with 
msw fit. For vacuum oscillations, the fit gives 5u ~ 10 -10 eV 2 and sin 2 2a; = 0.75. In the 
present analysis, 5 12 < 10~ 4 eV 2 , which is consistent with both msw and vacuum solutions. 
For u>, therefore, we choose two possible values, viz., u> small and uj large. These two typical 
choices cover the extreme ends of the possible effects of mixing in supernova neutrinos. 

If u is small (corresponding to the MSW solution to the solar neutrino problem), then 
the antineutrino survival probability becomes, P ee = cos 2 0cos 2 u; — > cos 2 — > 1, since <fi 
is small. This, in effect, is similar to the no-mixing solution. The large angle solution 
allows for a near-maximal mixing of uj ~ 45°; in this case, the survival probability becomes, 
P ee = cos 2 4> cos 2 uj — > 1 /2 and this corresponds to maximal mixing in the antineutrino 
sector. Therefore we have, 



1 
2 



F- ~ - f F- + F c , 



2F W ~ i (3F° + Fl) . (30) 

In any case, we have the result that P ee ^ 0.5 for any choice of uj when <p is small. Hence, 
typically, the antineutrino fluxes that reach the earth are combinations of V e and V x fluxes. 
Again, since the average energies of V e and V x are 15 and 24 MeV respectively, this results in 
an enhanced V e event rate and a reduced V x rate at the detector. It is important to note that 
these flux mixings are energy independent. For example, the energy spectrum of a given 
neutrino flavour produced in the supernova is not altered during propagation; however, its 
flavour content at the detector will depend on the extent of mixing. We shall now probe the 
quantitative effects of these mixings on the observed event rates. 

We first consider the case where there is no neutrino mixing. The largest contribution 
comes from the V e p interaction, which has a cross-section proportional to the square of the 
antineutrino energy. This is shown in the left-hand part of Fig. 1, where the number of 
events, N(k), in the k th bin is plotted against the central values of the recoil electron energy 
in that bin. In comparison, the V e e contribution is negligibly small. However, this is not 
the case with the V e O contribution, which though small, may be measurable at, say, the 
large (32 KTon) SuperK detector. It can be seen, though, that the total rate is saturated by 
the proton interaction. Neutrino mixing causes an increase in the high energy event rates, 
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as can be seen from the right hand side of Fig. 1. Here, we have used typical values for u 
(= 45°) and (= 10°); this choice of u maximises mixing effects. The other (small-angle) 
solution for u is similar to the no-oscillation scenario shown on the left hand side of the 
figure. Note that our three-flavour analysis precludes the choice u, <fi = 0. For example, the 
result that v e starts out as a pure v% mass eigenstate in the stellar core will not hold if <fi = 0. 

The corresponding results for v e events are shown in Fig. 2. While the no-mixing contri- 
butions are negligible, it is seen that there is a more than 10-fold increase in the event rate 
due to scattering off oxygen. The low no-mixing rate was because the average u e energy is 
less than the threshold energy required for this reaction to proceed. Mixing opens up this 
channel since there are now many more u e , originating as v x in the star, which are more 
energetic. Since the backward peak in the v e O cross-section is more pronounced for flux 
distributions at higher temperatures |25|, it may be possible to separate these events from 
the bulk of the anti-electron neutrino events at the detector. 

Finally, we see from Fig. 3 that there is a low-energy enhancement of the u x and V x 
rates upon mixing. Their contribution, however, is still small, of the order of the v e e elastic 
scattering events. For comparison, all the contributions, with and without mixing, are 
shown in Fig. 4. It is clear that the proton absorption events are the largest, independent 
of mixing. However, it is the v e O events which are most sensitive to the amount of mixing, 
and are likely to be most important in furthering our understanding of neutrino oscillations 
in vacuum and matter. 

We recall that the proton absorption events are isotropic, the scattering off oxygen is 
backward enhanced, while the elastic scattering on electrons is mostly forward peaked. In 
fact, even the elastic v e and v x (including antineutrinos) events may be separated based on 
angular resolution |2(J. Hence, it is likely that a nearby supernova explosion (at a distance 
of about 10 KPc, say) can yield information independently on the various neutrino flavours, 
z/ e , V e and v x . We have therefore shown the total contribution from each of these flavours, 
with and without neutrino mixing, in Fig. 5. For a 32 KTon detector such as SuperK, this 
translates to a total event rate of 12,235 events with mixing as opposed to 9,871 events 
without mixing, a 25% increase, with individual channels contributing as shown in Table 1. 



Note that the thermal flux distribution, while agreeing with the numerical model of Ref. (To 
at lower energies, overestimates the flux at larger energies. Hence, the number of events at 
high energies may be overestimated in this model. However, we emphasise that the relative 
increase, with and without oscillation, remains the same. 

We now briefly discuss these results in relation to the supernova SN1987A. Recall that 
the supernova, which was 55 KPc away, was detected by KII, which was a 2.14 Kton (fiducial 
volume) water Cerenkov detector; the corresponding results for this can therefore be obtained 
from our analysis by multiplying the results by a factor of 2.14/30.25. However, since KII 



mostly detected low energy events, we have now included the detector efficiency |17[ in our 
analysis. The event rate is determined entirely by the V e -p events. Since KII measured 
the time dependence of the spectrum, we have shown our results for the event rate as a 
function of time in Fig. 6. The figure on the left hand side of Fig. 6 shows the event rates 
for energies from 8 < E e (MeV) < 30 MeV, which is the energy range in which KII made 
observations; the dotted curve indicates the contribution in the absence of mixing. It is seen 
that mixing marginally decreases the event rate while it almost doubles the rate for events 
with energies E e > 30 MeV as can be seen from the figure on the right. This separation has 
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been done since the thermal Fermi-Dirac distribution that we use overestimates the flux at 
larger (neutrino) energies due to a very large high-energy tail compared to the corresponding 
numerical model |15|], while agreeing quite well at lower energies, E Pe ~ E e < 30 MeV. This 
is particularly true for the event rate at early times; hence our high energy predictions may 
be overestimated by a factor of 4 or more. Note, however, that even if the absolute spectrum 
is overestimated by the model we have used, the results we had shown earlier contrast the 
relative differences with and without mixing and still hold. Finally, the model also predicts 
that high energy events are most likely to occur at early times in the supernova explosion. 
This is not inconsistent with the observations of KII |J. 

To summarise, a great deal of the physics of neutrino mixing and the effects of dense 
matter in neutrino propagation may be tested from neutrinos emitted during supernova 
explosions. Specifically, from the results summarised in Table 1, we may conclude, 

1. The observed T7 e p events are the largest in number as well as least sensitive to the 
mixing parameters. Hence they provide a direct test of the supernova models. Since 
the angular distribution of these events is isotropic, they may be used to set the overall 
normalisation. 

2. All the interactions involving electrons as targets are peaked in the forward direction 
(in fact, for E u > 8 MeV, more than 90% of them lie in a 10° cone with respect to 
the supernova direction). In the absence of any mixing, there will also be a few events 
in the backward direction due to CC scattering on oxygen targets. As indicated in 
Table 1, the forward-backward asymmetry in the event distribution will be clearly 
marked. 

3. The main effect of mixing is then to produce a dramatic increase in the events involving 
oxygen targets. As remarked earlier, this will show up as a marked increase in the 
number of events in the backward direction with respect to the forward peaked events. 
The actual increase, however, will depend sensitively on a combination of both the 
mixing parameters as well as the supernova model. 

We have limited our analysis in this paper to a model with three active neutrino generations 
and possible mixings among them. This allows us to incorporate constraints arising from 
solar and atmospheric neutrino problems. This however leaves out a new set of constraints 



which may emerge from the results of the LSND experiments |31| . The LSND results cannot 
be accommodated within the three generation formalism, with the parameter ranges used 
in our analysis. One may therefore require a sterile neutrino with a new mass scale, leading 
to yet another mass squared difference, S ~ 0.3-1 eV 2 [^. In the context of supernova 
neutrinos, this opens up yet another channel for T to oscillate and may therefore reduce 
the dramatic enhancement one sees in the v e events. 

We thank Kamales Kar for many useful discussions and a critical reading of the paper. 
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APPENDIX A 



The time evolution of the neutrino mass eigenstates in vacuum is given by, 

Ui(t) = exp[-iEit]vi(0) ;i = 1,2,3. 
Assuming the neutrino masses to be small, in the extreme relativistic limit, we have, 

2p 

where //j, (i = 1,2,3), denotes the neutrino masses. In the presence of matter, neutrinos 
interact with electrons, protons and neutrons in matter. While u e (v e ) interact both via CC 
and NC interactions, u x (v x ) scatter via NC interactions alone. Note that interaction with 
matter is diagonal in the flavour basis but not in the mass basis. As a result, the dispersion 
relation in matter is given by, 

Ei ~ P + -± , 

2p 

where mj are now eigenvalues of the (mass) 2 matrix given by 

M 2 = M 2 + M- 2 + . 

m v 1 mt 



Here the mass-squared matrix in vacuum, 
given by 



M 2 , is defined in Eqs. (3) and (g) and M int is 



M 2 



±V2G F p [- {(1 - 4 sin 2 6 w )(N e - N p ) + N n ] I + 2A e M A 



Here N e , N p and N n denote the number densities of electrons, protons and neutrons in 
matter, Ma is the matrix defined in Eq. (|7|) of Sec. 2, and 6w is the Weinberg angle. Note 
that (1—4 sin 2 9w) is close to zero and that the matrix, M 2 nt , is expressed in the flavour basis. 
The upper sign corresponds to neutrinos (z/ e , and u T ) while the lower sign corresponds to 
antineutrinos (y e , V T ). 

We will now compute the eigenvalues, m 2 , in matter, ignoring the terms proportional to 
the unit matrix, I, since we will only be interested in differences between the eigenvalues. 
Consider the eigenvalues of the matrix defined by 



AM4 = A(r)M A + 5 31 M 3l + 6 2l M 21 



where 5ij = fi 2 — fi 2 and M 3 i, M%\ are matrices defined in Eq. (|j) and A(r), the matter 
dependent (and hence distance dependent) term, A(r) = ±.2\^2GpN e E for v e and W e re- 
spectively, varies linearly with the matter density. Since neutrinos are produced in the high 
dense region of the supernova, |A(core)| ~ 2 x 10 7 E eV 2 where E is the neutrino energy in 
MeV. Thus |A(core)| > <J 31 > 5 21 . 

We now compute the eigenvalues perturbatively. Since 8%i and are different from each 
other, the resonances, if they occur, are well-separated. We therefore diagonalise the first 
two terms in AM^ and treat the third term as a perturbation. The eigenvalues are then 
given by, 
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5m\ 



A + 5; 



31 



+ 2 



(A - 5 31 cos 20) 2 + (5 31 sin 2(j 



1/2 



+ 6 2 i cos 2 (0 - m ) sin 2 



5m 2 , = S 2 i cos 2 ; 



5m 2 



A + 5: 



31 



1/2 



2 2 

where the matter mixing angles are given by 



(A — 631 cos 20) + (5 3 i sin 20) + 5 2 i sin (0 — m ) sin ui 



(31) 
(32) 
(33) 



and 



tan 20 r , 



tan ax, 



5 3 i sin 20 
S 3 i cos 20 — A 



To the leading order, the mixing matrix in matter, C/ m (0 m ,u; m ), is given by, 



As, 
1 -A 

S 4>m AC0 m C0 m 







which is unitary, up to 0((^2i/^3i) 2 ) • Here c and s stand for cos and sin respectively, 
for example, c<p m denotes cos0 m ; A = {5 2 i/ hijs^^s^. Terms of order 0(5 2 \/A) are 
neglected in U m . 

Uptil now, the only approximation that has been used is the hierarchy, A ^> 5 3i ^> 5 2 i- 
Using this hierarchy and the value of A (core) given earlier, we find, 



(f>r, 



7r/2 for neutrinos ; 
for antineutrinos. 



This result, when combined with the definition of the matter mixing matrix, U m , leads to 
the fact that z/ e 's are produced almost entirely in the ^3) mass eigenstate whereas V e are 
produced almost entirely in the mass eigenstate, \Vi), in the dense stellar core. 

If the propagation is adiabatic, this also implies that the averaged survival probability 
of electron neutrinos is given by 

P ee =\(u e (t)\u 3 (0))\ 2 = sl, 

and that of the anti-electron neutrino is given by 

P ee =|(77 e (t)|77 1 (0))| 2 =4c 2 , 

where t denotes the time of detection of the neutrino on earth. 

In Appendix B we show that the propagation is indeed adiabatic and the expressions 
given above provide a reasonably accurate description of the matter effects in the stellar 
interior. 

A few remarks about the eigenvalues, Smf, are in order. Note that the eigenvalues 
themselves are always positive definite for electron type neutrinos, whereas for muon type 
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neutrinos or electron-antineutrinos, this is not always the case since A is large and negative. 
The complete dispersion relation for the energy eigenvalues is given by 

Ei ~ p + ±- f/xf =F V2G F p {(1-4 sin 2 6 w ){N e - N p ) + N n \ + bm\ 
ip 

The effect of the CC interactions with matter gives rise to the Smf term. The NC term, 
common to all flavours, is now included here. The second term in the expression for E{ is 
typically of the order of tens of eV in the stellar core. Therefore, for energies E ~ p ~ 
few MeV, the second term is small and may be neglected except when computing matter 
mixing angles. However, for neutrinos having energies of the order of tens of eV (but 
still with p ^> Hi), the two terms compete. Since the sign of the second term changes 
depending on whether the particle scattering is a neutrino or an antineutrino, one may 
expect interesting phenomena when E^ becomes negative. This may lead to the trapping 
of low-energy neutrinos. This phenomenon is unique to neutrinos produced in supernova 
explosions. While it is of little relevance to the detection of neutrinos on earth, it may have 
interesting astrophysical consequences. The dynamics of such neutrinos is under further 
investigation. 



APPENDIX B 



The electron neutrino v e is produced in the core of the supernova in the mass eigenstate 
1 1/3) with a negligible admixture of the other two states. As the produced \v$) propagates 
outwards, it passes through variable density matter (since the density inside the stellar 
core decreases monotonically outwards). Such a propagation may in general induce the 
presence of other mass eigenstates (since they are no longer eigenstates of the Hamiltonian) 
by the time the neutrino exits the star and reaches the detector. The Landau-Zener "jump 
probability" or level transition probability P§[ is maximal at the resonances and is given by 



LZ 



cxp 



7T 

7 F 

2 ' 



where F is a factor which depends on the density profile and 7 is the non-adiabaticity 
parameter. 

Since v e is produced in the mass eigenstate |z/ 3 ), we first consider the cross-over between 
1 1/3) and 1 1/2). Then, 

5 sin 2 2<p 



7 



IE cos! 



1 dN P 



N P dx 



where <fi is the relevant mixing angle for the upper resonance and 5 = (£31 + £32 + 
<5i2 cos2u;)/2 ~ <5 31 independent of cos2u; |L6|], because of the assumption, 5u S31 ~ 632- 



The suffix '0' indicates that the derivative in the density, N e , is to be evaluated at resonance, 
when the eigenstates \u 2 ) and |i/ 3 ) are the closest. 

The density profile in the core may be assumed to be of the form [119, 



p(r) 



C 

~3 ' 



r 



15 



where 1 < C/10 31 gms < 15. (This assumption is not crucial but is sufficient for our analysis). 
As a result, the non-adiabaticity parameter evaluates to 

R0831 ( sin 2 20 \ 

7 



6E \ cos 20 I 

where i?o is the radius at which the higher resonance occurs, i.e., when A(r) ~ £31 cos 20. 
Using the explicit expression for A(r), the resonant density is given by 

, . N 6.6 x 10 5 
Po(gm/cc) = 

where £31 is in eV 2 and e = E/(10 MeV). Here Y e = Z/A ~ 0.5 is the electron fraction in 
the matter. For £31 ~ 10~ 3 eV 2 , as preferred by solar and atmospheric neutrino data, and 
a typical detected neutrino energy of 10 MeV, the resonant density is po = 132O(cos20) 
gm/cc. This determines the resonant radius, Rq for a given value of 0. If the angle is 
indeed small, as indicated by chooz |^7|, this implies 




/ eC \ 1/3 

R = [ ) km , 

° V1320 xl0 15 / 

which evaluates to 20-50 thousand kilometres for E = 10 MeV (e = 1) 
The non-adiabaticity parameter is then given by (for Ro in km), 



7 



5076i?o /sin 2 : 



cos 20 



Since Rq is large, it is clear that 7 is large unless is very small. Furthermore, for small 
values of 0, F ~ 1. In fact, for sin0 > 10~ 2 , we find 7 ^> 1 so that the Landau-Zener 
probability is vanishingly small: P^z < 10 -2 . Our three-flavour analysis of the neutrino 
mixing problem in any case precludes the choice of = 0. Therefore, for all practical 
purposes, we assume > 10~ 2 and hence consider the neutrino propagation in matter to be 
purely adiabatic. This implies that the \v e ) which is produced as \v%) essentially remains in 
this mass eigenstate until it reaches the detector. For I7 e , which are produced mainly in the 
\V\) mass eigenstate, there is no resonance condition to be satisfied (the sign of A changes 
from neutrino to antineutrino) and hence the propagation is always adiabatic. 
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TABLES 



Detector 


v e e 


v e (J 


v e e 


VeO 


v e V 














E e >8 


MeV 






1 KTon (no osc) 


2.3 


1.0 


0.8 


3.8 


272.0 


1.3 


1.0 


1 Klon (max oscj 


4.2 


23.8 


0.9 


8.1 


323.2 


1.0 


1.0 


SuperK (no osc) 


72.4 


30.8 


25.1 


123.0 


8702.9 


41.1 


33.4 


SuperK (max osc) 


134.6 


761.0 


30.0 


260.6 


10343.9 


31.5 


31.7 










E e > 30 MeV 






1 KTon (no osc) 


0.0 


0.2 


0.0 


1.8 


76.6 


0.2 


0.2 


1 KTon (max osc) 


0.7 


18.4 


0.1 


5.4 


138.6 


0.1 


0.1 


SuperK (no osc) 


1.5 


5.9 


1.0 


56.2 


2450.0 


6.4 


4.9 


SuperK (max osc) 


22.2 


587.6 


2.2 


172.8 


4436.5 


3.5 


4.1 


TABLE I. Event rates with and without oscillation for 


a supernova explosion 


at a distance of 


10 KPc. The high 


energy events, 


which 


are given 


separately, show the 


enormous 


enhancement in 



the p and O channels, with oscillation. Results are shown for a 1 kTon water Cerenkov as well as 
for the 32 kTon (fiducial volume) SuperKamiokande detectors. 
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FIG. 1. The number of events in bins of electron energy of 1 MeV each, due to V e interac- 
tions, are shown as a function of the electron energy, with and without mixing. The long-dashed, 
dashed, and dotted lines correspond to interactions with p, O and e respectively, in the detector. 
The dot-dashed line indicates the effect of inclusion of detector efficiency and resolution on the 
interaction with p. See text for more details. The solid line denotes the total contribution to the 
event rate from V e . 
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FIG. 2. The number of events in bins of 1 MeV each, due to u e interactions, are shown as 
a function of the electron energy, with and without mixing. The dashed and dotted lines corre- 
spond to interactions with O and e respectively, in the detector. The solid line denotes the total 
contribution to the event rate from u e . 
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FIG. 3. The number of events in bins of 1 MeV each, with and without mixing, due to f MjT e 
and Vfi )T e elastic scattering in the detector, are shown as a function of the electron energy, as 
dotted and dashed lines respectively. The solid line denotes the total contribution to the event rate 
from all these channels, that is from v T , V^, and V T . 




FIG. 4. A comparison of the number of events in bins of 1 MeV each, due to various processes, is 
shown as a function of the electron energy, with and without mixing. The line types indicate events 
from the processes V e p (solid), V e e (dotted), V e O (dashed), u e e (long-dashed), v e O (dot-dashed), 
and v x e (dot-long dashed) processes respectively. The subscript x denotes the NC contribution 
from u T: and their antiparticles. 
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FIG. 5. The total number of events (summed over all processes) in bins of 1 MeV each, due to 
X7 e , z/ e , z/„ T and Vu T interactions, are shown as a function of the electron energy. The solid and 
dashed lines denote the event rates without and with (maximal effect due to) mixing. 
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FIG. 6. The time-dependent neutrino spectrum (due to V e p scattering) in bins of 0.5 sec is 
shown as a function of the time of detection, in comparison with the events observed at the 
Kamiokande II detector from the supernova SN1987a. The dashed and solid lines correspond to 
the number of events without and with (maximal effect due to) mixing. The low and high energy 
components of the signal are separately shown. 
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